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DIS EVENT-SHAPE RESUMMATIONS AND SPIN-OFFS 3 
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We present results from a recently completed project to calculate next-to-leading logarithmic 
resummed distributions for a variety of event shapes in the 1 + 1-jet limit of DIS. This allows 
fits for the strong coupling and for non-perturbative effects using the large amount of data on 
these observables from HERA. Spin-offs include the discovery of a new class of logs for certain 
final state observables (non-global observables); a program that allows a speed-up by an order 
of magnitude of certain fixed-order calculations in DIS with DISENT or DISASTER++; and 
the development of state-of-the-art PDF evolution code. 



1 Introduction 

Event-shapes are observables sensitive to the flow of energy and momentum in hadronic final 
states. They have been extensively studied in e + e~ collisions, for example for the measurement 
of the strong coupling, tests of QCD through fits for the colour factors and the study of novel 
approaches to hadronizationS Typically the most discriminatory studies make use of event-shape 
distributions, which are compared to next-to-leading perturbative predictions that are resummed 
in the 2-jet limit. 

Recently the HERA experiments have also started considering event shapes, defined in the 
current hemisphere of the Breit frame. In particular, distributions have been measured by Hl@ 
and while only mean values have been so far studied by ZEUsi it is our understanding that 
they intend to extend their studies to distributions. 

Resummed predictions for e + e~ event shapes have existed in the literature since the early 
nineties^ but until recently no such calculations were available for DIS observables (an exception 
is jet rates!!). Because of the strong similarity between a hemisphere of an e + e~ event and the 
current hemisphere of the DIS Breit frame, it is natural to assume that the extension from 
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e + e~ to DIS will be fairly straightforward. It turns out not to be so, for both conceptual and 
technical reasons. In what follows we outline some of the issues that arise, and then present 
some preliminary comparisons to data. 

2 Resummation issues in DIS 

The most obvious new issue to arise in DIS compared to e + e~ is that of collinear factorization. 
Despite the fact that the observables are all defined in the current-hemisphere (Tic) of the Breit 
frame, owing to details of the kinematics those defined with respect to the photon axis are 
sensitive to emissions in the remnant hemisphere (Tin) through recoil effects. Requiring the 
event-shape to have a value close to that of the 1 + 1-jet limit, one therefore forbids emissions 
in the whole of the phase space (Tic an d Tin). However collinear factorization at the scale Q 2 is 
conditional on there being no restrictions on emissions in the remnant hemisphere. Requiring the 
event shape to have a value less than some V, which translates to a limit on the largest possible 
transverse momentum of collinear emissions in Tin, &t ^$ V n Q 2 (n is observable-dependent), has 
the consequence that collinear factorization can only be recovered if parton distributions are 
evaluated at a scale of the order of V n Q 2 . This is actually a familiar result from calculations of 
the pt distribution of Drell-Yan pairsi and has also been observed in more complicated multi-jet 
event shapes! 




Figure 1: Contributions relevant in the calculation of non-global terms; the triple lines indicate incoming partons. 

A second issue is that of non-global logarithms, which arise in observables sensitive only to 
emissions in a restricted portion of phase space (e.g. Tic) such as the jet mass, and additionally 
in observables whose sensitivity to emissions is discontinuous across one or more boundaries 
in phase space (an example El is the thrust t z e)- An erroneous assumption that has widely 
been made in the literatureH'li30t3 is that (to single-logarithmic accuracy) in order to suppress 
radiation into Tic (say), it suffices to suppress primary radiation from the various hard 'legs' 
into that hemisphere. 

While at leading order in a s In V this is correct, starting from second order, configurations 
such as those in fig. |l] become relevant. El The crosses indicate emissions which must be forbidden 
in order for the observable to have a small value. The grey emissions are those that do not 
directly affect the value of the observable. The left-hand picture represents the configuration 
relevant at second order: a soft emission (1) in Tin, which does not contribute to the observable, 
radiates an even softer emission (2) into Tic, which does contribute to the observable. The strong 
ordering in energies Q > ]?i » i?2 leading to one power of InV for each power of a s . While 
this term is calculable analytically, at all orders one needs to forbid coherent radiation into Tic 
from arbitrarily complicated ensembles of large-angle energy-ordered gluons in Tin (right-hand 
picture). This is complicated both from the point of view of the colour structure and of the 
geometry. The former can be dealt with approximately in the large Nq approximation, while 
the latter can so far only be treated numerically. Some insight into the dynamics associated 
with these non-global logs was obtained in the context of a more general study of energy flow 
distributions where one finds that in the limit of large In 1/V not only is radiation into Tic 
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Figure 2: left, a comparison between resummed predictions and HI results; right, results of fits for a s and ao- 



forbidden, but radiation at intermediate energy scales is also forbidden in a neighbouring 'buffer' 
region of Ti-ji- The size (in rapidity) of this buffer region increases with lnl/ and the overall 
suppression factor coming from non-global logs seems, at least in part, to be associated with the 
suppression of primary radiation into the buffer. 



3 Technical issues 



When implementing the resummations as computer programs to allow comparisons to data a 
number of technical issues arose. When this project started there existed only two subtraction- 
based programs for NLO calculations in DIS, namely DISENT0 and DISASTER++0, which 
were known to disagree for certain observables.HZI Comparisons with the expansion of the re- 
summed results made it possible to identify DISASTER++ as the one giving correct predictions^ 

Unfortunately, of the two programs, DISASTER++ is an order of magnitude slower, and we 
would have needed over a year's computing time to obtain the fixed-order predictions (needed to 
describe the observable in the region outside the 1 + 1 jet limit) to sufficient accuracy. However, 
the traditional approach to such calculations involves considerable duplication of effort: one 
needs distributions at several x and Q 2 values and typically one uses separate events for each 
x and Q 2 value. But the matrix elements (modulo their ysj -dependence) and the calculation 
of the event-shape are both independent of x and Q 2 , so each event in an NLO Monte-Carlo 
program can be 'reapplied' to several x and Q 2 values. We have written a program DISPATCH, || 
which acts as a wrapper to DISENT and DISASTER++ so as to automate such a procedure. 

Another spin-off from this project is the development of a high-precision PDF evolution 
code^l which has been used in collaboration with VogtH to produce reference NNLL evolutions 
to an accuracy of 1 part in 10 5 . 



4 Comparison to data 
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land plot of fig. gl shows a comparison between our matched resummed distribu- 
1 and the HI dataH for the C-parameter, Ce- Non-perturbative contributions have 
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been included using 1/Q corrections, whose size have been hypothesized Ell to be governed by 
a universal parameter olq. We have fitted for both a s and otQ, using only the red points of the 
distributions. The results (1-<t contours) for Ce and a number of other DIS observables are 
shown in the right-hand plot (green curves) and compared to e + e~ resulte for mean values (red 
curves), with the jet masses measured in so-called massless schemes. B3 The agreement both 
within DIS and across experiments is strong confirmation of the universality hypothesis for ao- 
There remain some observables in DIS where the agreement is less good, notably those 
measured with respect to the photon axis (t z e, B z e), and the situation worsens if one includes 
lower Q (also lower x) data. The detailed origin of the problem remains to be understood, 
though it may well be associated with higher-order corrections being relatively larger at lower 
Q values. In this context, higher precision data (and higher resolution data as well), especially 
at larger Q values, will be interesting as it will make it possible to pin down any systematic Q 
dependence over and above that expected from the theoretical predictions used so far, perhaps 
for example from shape functionsCJ 
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